Explicit illustration of non-abelian fusion rules in a small spin lattice 
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We exactly solve a four-site spin model with site-dependent Kitaev's coupling in a tetrahedron by 
means of an analytical diagonalization. The non-abelian fusion rules of eigen vortex excitations in 
this small lattice model are explicitly illustrated in real space by using Pauli matrices. Comparing 
with solutions of Kitaev models on large lattices, our solution gives an intuitional picture using 
real space spin configurations to directly express zero modes of Majorana fermions, non-abelian 
vortices and non-abelian fusion rules. We generalize the single tetrahedron model to a chain model 
of tetrahedrons on a torus and find the non-abelian vortices become well-defined non-abelian anyons. 
We believe these manifest results are very helpful to demonstrate the nonabelian anyon in laboratory. 
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The spin lattice models of Kitaev-type have attracted 
many research interests because of the abelian and non- 
abelian anyons in these exactly soluble two-dimensional 
models 0, Q , which are of the potential application in the 
topological quantum memory and fault-tolerant topolog- 
ical quantum computation 0]. 

The abelian anyons can be explicitly shown in Kitaev 
toric code model [l[ or Levin- Wen model Fusion 
rules and braiding matrix can be easily illustrated in real 
space. This has simulated many attempts to design and 
process experiments to demonstrate these abelian anyons 
in laboratory 

In solving Kitaev honeycomb model [l| and its ram- 
ifications [a S H]) a key technique is the usage of the 
Majorana fermion representation of the spin- 1/2 oper- 
ators. The systems then are transferred into bilinear 
fermion systems and the ground state sector can be di- 
agonalized in the momentum space. The shortcoming to 
use the momentum space is that the ground state and 
the elementary excitations are hard to be expressed by 
the original spin operators, i.e., Pauli matrices. On the 
other hand, the sectors with vortex excitations can only 
be treated numerically. Then the nonabelian fusion rules 
and statistics may not be directly shown in Pauli matri- 
ces' language. Lahtinen et al have derived the nonabelian 
fusions through the spectrum analysis 0]. However, they 
are still not directly related to Pauli matrices. A real 
space form of the ground state of the Kitaev honeycomb 
model in the insulating phase with abelian anyon was 
studied [loj | . We tried to present the non-abelian fusion 
rules for high energy excitations in this Kitaev model 11]. 
There was an attempt to use toric code abelian anyons to 
superpose the Ising non-abelian anyons without involving 
a Hamiltonianflil] . 

In this paper, we solve a spin model with Kitaev's cou- 
pling in a small lattice, i.e., a tetrahedron (Fig[T|). We 
perform Majorana fermions and their zero modes, non- 
abelian vortices and their fusion rules in real space by 
means of Pauli matrices. Because the system is finite and 
everything can be deduced in an elemental way, it will be 



very helpful to intuitionally understand these concepts 
which were used to be expressed in those deep mathe- 
matical language. Since the spin configurations of these 
excitations and the fusion rules are explicitly shown, the 
experimental techniques with cold atoms, finite photon 
graph states and nuclear magnetic resonance systems [f| 
are possible to be applied to demonstrate these quantum 
states and then to design quantum bits and gates for a 
topological quantum computer. 

The single tetrahedron model is too small for non- 
abelian anyons to be well-defined. We generalize it to 
be a chain model of tetrahedrons on a torus. This chain 
model is also exactly solvable and the non-abelian anyons 
can be well-defined. They are still of a simple form like 
the non-abelian vortex in the single tetrahedron model. 




FIG. 1: (a) A tetrahedron in whose points 1,2,3 and 4, the 
spins live. The four surfaces are labelled by A, B, C and D or 
124, 134, 234 and 123. (b) The top view of the tetrahedron. 
x,y,z are the links with different Kitaev couplings. 



The model and symmetries Kitaev model in a tetrahe- 
dron is given by H K = Y, x •/.-.,.,« + Y, y Jy,ij a iVj + 
J2z Jz,ij a i a j where x,y,z are the links shown in Fig. 
[IJb) and J a ,ij are link- and site-dependent coupling con- 
stants, a" are spin-1/2 operators obeying Pauli matrix 
algebra, e.g., of of = iof and (erf) 2 = (erf) 2 = (of) 2 = 1. 
The spin operators on the different sites commute, i.e., 
[crf,Oj] = for i ^ j. An intuitive imagination is the 
model may be easily solved if J a ,ij are not site-dependent 
as that in Kitaev model in an infinite lattice Q|. However, 
a direct check finds that, unlike Kitaev model in an infi- 
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nite lattice, this model can not be reduced to a bilinear 
fermion theory in such a coupling constant choice. This 
is because the tetrahedron is topologically equivalent to 
a sphere, Kitaev model defined on this compact space 
is very different from the model on infinite or periodic 
lattices. In this paper, we consider these coupling con- 
stants are site-dependent. We also include some three- 
and four-spin coupling terms as those in a generalized 
Kitaev model jf|. The model Hamiltonian we will study 
is given by 



The eigen vaules of the Hamiltonian matrix are 
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There are four conserved operators which live in triangu- 
lar plaquettes: 

P={P A = Pl24 = <jlo\ol, P B = Pl34 = 0?<T%*1, 

Pc = P234 - of Pd = P123 = <r\°2°l], (2) 

which arc mutual commutative. They obey [Pijk,H] = 
and PaPbPcPd = 1- Three-spin couplings break time 
reversal symmetry. The Hamiltonian is time-reversal in- 
variant if k = 0, . 

Bilinearization, diagonalization and states This Hamil- 
tonian can be transferred into a bilinear fermion Hamil- 
tonian. In the deduction process, we can illustrate some 
abstract concepts in an elementary way. For example, 
we can define eight Majorana fermions corresponding to 
four sites: 



4>i = 01,^2 = ofof,^ = 03°2 cr li ^4 = alala z 2 a\ (3) 
h = -of ,62 = -o-jfofj&a = -ofofo-f,^ = -a\alala\ 

They obey {ipi,i/)j} = {bi,bj} = 2<% and {ipi,bj} = 0. 
Remarkably, [H, bi\ = 0. This can be directly checked or 
be seen by writing the Hamiltonian in terms of Majorana 
fermions, 

H = -t{d\d b + d\d a ) - f i(dld a + dld b ) + A*d a d b + Adldl 

where d a — h(ipi + ^2) an d d b — ^(^3 + iip4,)- The 
parameter relations are given by fi = 2J X , t = J y + 
A, A = Ai+iA2 = Jy—X+in. An easy way to identify ([T|) 
and ([4]) is substituting ([3]) into ([4]) . One can check d a and 
db are conventional fermions, i.e., {d a ,d\} = {d b ,d b } = 
1 and {d a ,d a } = {d b ,d b } = {<4,d&} = {d a ,d\} = 
[q, 1 1 01 ] . This is a BCS p-wave pairing Hamiltonian in a 
finite system and can be diagonalized by rewriting Q as 
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E = --(^/WT7 + t),E t = _-( v /jAja+7?-t), 



= -(vW+A?-*),-^ = -(vW+? + t), (6) 
Diagonalizing the Hamiltonian, one has 

H = E 2 d\d b + E^J a +E Q d b d\+E t d a d\. (7) 



The generalized Bogoliubov fermions d ajb can be ob- 
tained in a standard way with d a>b — %\^d a fi + 

x 2 a ' b ^d bta + £3 b + x^^dl a and the coefficients x^ ,b 
are normalized eigen vector of the Hamiltonian ma- 
trix. The Bogoliubov fermion operators obey the stan- 
dard fermion commutation relations. A subspace of 
quantum states are {\G), d a \G), d\\G), d\d\\G)} where 
\G) = d a d b \Q) with d„|0) = 4|0) = 0. The vacuum 
|0) = d a d b \<f)) for a reference state \(j)), e.g., | T T T T ) - 
The eigen energies of this set of quantum states are 
{-VIAP + M 2 , -t, t, v/|A| 2 + M 2 }. When ^\A\ 2 + ^ > 
t, i.e., 4J 2 + k 2 - AJ y \ > 0, \G) is the ground 
state. Because [if, bi] = 0, each energy level is for- 
mally sixteen-fold degenerate, e.g., the ground states are 
{|G>,6i|G>,6i6 i |G> ) 6 i 6 i 6 fc |G>}. That is , h play a role of 
zero modes of Majorana fermions. However, there are 
only four independent, which, e.g., are 



{\G),c\\G),4\G) 

and Co 



\c\\G)} 



(8) 



where c[ = h(b% — ib^) and c\ = \(b 2 — ib^) with 
C\\G) = C2\G) = 0. The total Hilbert space is sixteen- 
dimensional as expected. Four degenerate states in a 
given energy level are distinguished by quantum number 
P = (Pa, Pb, Pc,Pd), which are shown in Tab. 1. 
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C 2 


J J 
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\G) 


(1,1,1,1) 


(-1,1,1-1) 


(1,-1,-1,1) 


(-1,-1,-1,-1) 


dl\G) 


(-1,-1,1,1) 


(1,-1,1-1) 


(-1,1-1,1) 


(1,1,-1,-1) 


dl\G) 


(-1,-1,1,1) 


(1,-1,1-1) 


(-1,1-1,1) 


(1,1,-1,-1) 


dt,di\\G) 


(1,1,1,1) 


(-1,1,1-1) 


(1,-1,-1,1) 


(-1,-1,-1,-1) 



Table 1: The eigen values of P of the quantum states. 

Due to the constraint PaPbPcPd = F there are eight 
different P which are carried by the states in the first two 
levels or in the last two levels as shwon in Tab. 1. 

Fusion rules: abelian and non-abelian We now go to 
illustrate the fusion rules of these eigen excitations. First, 
we define Majorana fermions corresponding to d a y. 



H = ±(4,dl,d a ,d b 



( -[i -t 
-t -n A 
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i>i = da + dl, lp2 = ~i(d a - d\), 



i>3 = d~b + d\, ip. 



-i(d b 



(9) 



They obey {V'iiV'j} = ^ij an( i {^iibj} = 0. There are 
(5) four sets of states which obey abelian fusion rules, as 
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those in Kitaev toric code model, 



^U 3 > ^fc-^, a^fc-^, (10) 



7 (1) 



for i — 1, 2, 3, 4. These operators are 



■ j 7 (2) (3) 
rJ, 7 ( 2 ) h ( 3 ) 



-#i; 
= #2; 



(ii) 



-&1^>1&303, <4 



-6361^1, ct 3 = fabi^i; 



O4 = -h^bi^A, Vf' = -&4&2^2, = --04^2^2- 

Acting on the ground state |G), they are eigen states of P 
and their eigen values can be read out from Tab. 1. They 
are also eigen states of the Hamiltonian and their eigen 
energies can be read out from the number of tpi in a given 
operator. Since each energy level is four-fold degenerate, 
we find that the linear combination of these degenerate 
states may obey non-abelian fusion rules. They can be 
thought as non-abelian vortices. For example, 

Va = —j=(b 2 ijj 2 - 6462^2), VB = —^{bitpi - 6 3 6iV>i), 



V2 X 



V2' 



They are in fact the superposition of those toric code 
abelian vortcies. (We will be back to this point when 
we study braiding of anyons.) Acting these operators 
on |G), they are also eigen states of H. The details of 
P and H's eigen values list on Tab. 2. The subscript 
indices are used because, e.g., va\G) has an eigen value 
Pa = — 1 and other three either are 1 or do not have a 
definite eigen value. Therefore, t\a can be thought as a 
vortex excitation on the surface A, and so on. 
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Table 2: The eigen values of P and //of rj. '*' means the 
vortex does not have a definite eigen value. 

Using the algebra of Pauli matrices or equivalcntly, 
the anti-commutation relations between the Majorana 
fermions, one may directly prove that these operators 
obey the following non-abelian fusion rules 

VAVA ~ 1 + bi, VAbi ~ TjAi bibi = 1, 

VbVb ~ I + 63, VbH^vb, hb 3 = I, (13) 
VcVc ?~1 + Z»i, Vch^Vc, Mi = 1, 
VdVd ~ 1 + b 2l VDh^VD, b 2 b 2 = l, 

which are standard non-abelian Ising fusion rules. Equa- 
tions (|13p are one of central results in this paper. Only 



when the zero modes of Majorana fermions exist, the 
non-abelian vortices are eigen excitations [l3j]. We hope 
this illustration can help condensed matter physicists 
have a direct impression to these elusive mathematical 
relations and understand them in an elementary way. 

Breaking of time reversal symmetry The Ising non- 
abelian fusion rules and time reversal symmetry are con- 
comitant. The three-spin coupling terms in Hamiltonian 
(fTj) explicitly break the time reversal symmetry. How- 
ever, our deduction of the non-abelian fusion rules does 
not rely on a non-zero n. They hold even n = 0. The 
only change is A2 = K = and the gap A = Ai. Actu- 
ally, when k = 0, the time reversal symmetry is sponta- 
neously broken. The ground states are four-fold degener- 
ate, which are given by ([8]). Since under the time rever- 
sion T, Ta a T- x = -a a , one has T\G) = \G), Tc\\G) = 
4\G), Tc\\G) = -c\\G), Tc\c\\G) = -c\c\\G). We see 
two sectors which have different eigen values of T, i.e., the 
time reversal symmetry is spontaneously broken. This is 
because of the geometric frustration of the tetrahedron. 
This spontaneous breaking of the time reversal symmetry 
first observed in Kitaev model on a triangle-honeycomb 



lattice and leads to a chiral spin liquid [14 1 



Non-abelian anyons A frequently quoted result is that 
the excitations obey non-abelian fusion rules like (fTB")) is 
equivalent to the braidings of these vortices 77 are non- 
abelian and then these vortices are called anyons with 
non-abelian statistics or non-abelian anyons [3j. How- 
ever, it is based on these anyons are well-defined and 
they are energetically degenerate. In this small system, 
to identify an n vortex as an anyon is reluctant because 
we see that, e.g., t\a is not an eigen state of Pb and Pc, 
which means this is not a particle-like isolated excitation. 
On the other hand, the vortices t/a.b and Vc.d are not 
in the same energy level if t ^ and braiding two vor- 
tices with different energy do not make sense in statistics. 
Therefore, for this small system, we only emphasize the 
non-abelian fusion rules of these vortex excitations but 
not call them non-abelian anyons. 




FIG. 2: left panel: A chain of tetrahedrons on a torus. Right 
panel: unwind the torus to a periodic lattice. The thick lines 
carry J x while the dash lines carry J y . 

To well define an anyon, we generalize the single tetra- 
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hedron model to a chain model of tetrahedrons on a torus 
(Fig. [2]). The model Hamiltonian is given by 



n 
i=l 



with n + 1 = 1. Only a half of triangular 
plaquette operators are conserved, which are P = 

{P2a-\ — P2a-l 1 2a-l 2 2a 1 = "fa-li °2a- 1 2 °2ai ' ^2a = 

f2o-ii2ai2a 2 = ^L-ia °2ai °fa 2 }• The Hamiltonian can 
be diagonalized as 



Missing of the complex phases e~™/ 8 in i?^ and e m l % 
in i?^ is because R = B) for the toric code anyons [l2| • 
A framing therefore is needed We do not intend to 
propose a framing and ancillary qubits to implement the 
non-trivial chirality but refer to Wootton et al [121] . 

In conclusions, we presented a simple model in which 
there are a set of vortices obeying non-abelian fusion rules 
which were explicitly illustrated in an elementary way 
without using deep mathematical tools. Finally, we gen- 
eralized the single tetrahedron model to a chain model of 
tetrahedrons on a torus and showed that the non-abelian 
vortices defined in the single tetrahedron model become 
well-defined non-abelian anyons. 

This work was supported in part by the national natu- 
ral science foundation of China, the national program for 
basic research of MOST of China and a fund from CAS. 



r chain = -t^2(dld i+1 + d\ +1 di) -/j^dt^ 

i=l i=l 

n 



i=l 



where di — 
are given by 



iipi 2 ) an d the Majorna fermions 
4 II, i°j 



rr^ rr z FT rr z rr z h 

"i2 <J h ilj<i u ji u ji' u h ~ llj<i"jx"j 

~ a i2 Cr hi[j<i cr ji cr ji> wnere Majorana fermions bi com- 
mute with -ffchain- Therefore, bibj • • • |) are degenerate 
states if |) is an eigen state. Since [P, i[>i] = for all i 
but not all bi, the eigen values of P for the vortex states 
T]i = -^(biipp — bjbiipp), j i are determined by bi,j. 

Here ip p are the linear combination of ipi after diagonal- 
izing -ffchain as those in the single tetrahedron model. If 
j is far from z, r\i has only one minus P near the i-th site. 
Then, rji is a well-defined single vortex operator and can 
be thought as a non-abelian anyon due to 77^?^ = 1 + bj. 
There are n such anyons which are degenerate. Each 
anyon brings a zero mode bi of Majorana fermion in its 
center 



13] 



Non-abelian braiding matrices Rewriting rji = 

rrii) with e; = ibiip p and = ibjbiipp, one has abelian 
fusuion rules em ~ b, mb ~ e, eb ~ m and e 2 = m 
b 2 = 1. e and m are toric code mutual abelian anyons 
with degenerate energy. Thus, the Ising anyon is in fact 
the superposition of the toric code abelian anyons. This 
result has been recognized in refs. Il|,|l2j but without in- 
volving the Hamiltonian. Therefore, the toric code braid- 
ing matrices determine the non-abelian braiding matrices 
of r\i [I2I ] , which are Ising- like braiding matrices 

< = -l, {RTf = 1, «') 2 = -l, {RT? = -i- 
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